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ABSTRACT 

Context. The physical structure of hot molecular cores, where forming massive stars have heated up dense dust and gas, but have not 
yet ionized the molecules, poses a prominent challenge in the research of high-mass star formation and astrochemistry. 
Aims. We aim at constraining the spatial distribution of density, temperature, velocity field, and chemical abundances in the hot 
molecular core GlO.47+0.03. 

Methods. With the Submillimeter Array (SMA), we obtained high spatial and spectral resolution of a multitude of molecular lines 
at diff'erent frequencies, including at 690 GHz. At 345 GHz, our beam size is 0.3'', corresponding to 3000 AU. We analyze the data 
using the three-dimensional dust and line radiative transfer code RADMC-3D for vibrationally excited HCN, and myXCLASS for 
line identification. 

Results. We find hundreds of molecular lines from complex molecules and high excitations. Even vibrationally excited HC^^N at 690 
GHz is detected. The HCN abundance at high temperatures is very high, on the order of 10~^ relative to H2. Absorption against the 
dust continuum occurs in twelve transitions, whose shape implies an outflow along the line-of- sight. Outside the continuum peak, 
the line shapes are indicative of infall. Dust continuum and molecular line emission are resolved at 345/355 GHz, revealing central 
flattening and rapid radial falloff" of the density outwards of 10"^ AU, best reproduced by a Plummer radial profile of the density. No 
fragmentation is detected, but modeling of the line shapes of vibrationally excited HCN suggests the density to be clumpy. 
Conclusions. We conclude that GlO.47+0.03 is characterized by beginning of feedback from massive stars, while infall is ongoing. 
Large gas masses (hundreds of M©) are heated to high temperatures above 300 K, aided by diff'usion of radiation in a high-column- 
density environment. The increased thermal, radiative, turbulent, and wind-driven pressure drives expansion in the central region and 
is likely responsible for the central flattening of the density. 
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1. Introduction 

Massive stars and star clusters are bom deeply embedded in 
molecular clouds (for a review, see [Zinnecker & Yorke 2007 ). 
When cores have sufficiently contracted to form massive stars, 
the dense dust and gas is heated by these stars. The ice mantles 
around dust grains evaporate, and a rich plethora of molecular 
lines can be observed , along with strong dust emission ( Kurtz| 
|et al.|[2000[ [Cesa roni 2005 ). Often, these hot molecular cores 
are associated with ultracompact or hypercompact Hn regions, 
which are ionized by newly formed massive stars ( |Hoare et al. 
|2Q07| ). The physical and chemical structure of hot molecular 
cores is of great importance for the study of high-mass star for- 
mation and of astrochemistry. Investigations of the structure are 
hampered, though, by the compactness of the sources, by their 
scarcity and large distance, and by the large foreground column 
density, which only long- wavelength radiation can pass. 

What is needed, hence, is high angular resolution at 
(sub)millimeter wavelengths, combined with high spectral res- 
olution of many molecular lines, which contain all the informa- 
tion about chemistry, velocity field, and temperature. With the 
Submillimeter Array (SMA) in Hawaii, we observed the mas- 
sive hot molecular core G10.47-h0.03 at around 200, 345, and 
690 GHz, yielding a best resolution of 0.3'' and 350 identified 
molecular lines. Due to its large dust mass and high tempera- 



tures, this hot core is among the strongest sources at submillime- 
ter wavelengths, although it is located at a distance of 10.6 kpc 
( Pandian et al. 2008 ). It h as an estimated luminosity of 7 x 10^ 
Lq ( [Cesaroni et al.|2QlQ| ) and displays excepti onally many lines 



Wyrowski et al. 



from highly excited molecules (e.g. HC3N, 
[T9991 ). 

2. Observations and data reduction 



The high-mass star-forming region G 10.47 -1-0.03 was observed 
with the Submillimeter Arra}j^(SMA) during five nights in dif- 
ferent array and receiver configurations, yielding six two-GHz- 
wide bands centered at around 201, 211, 345, 355, 681, and 691 
GHz. Table [T] summarizes the observations. We note that Jupiter 
was close (~ 10°) to the target source in 2008, so its moon 
Callisto could be used as gain calibrator at 690 GHz, where it 
is otherwise very difficult to find a suitable calibrator source. It 
was observed repeatedly for 2 minutes after spending 10 minutes 
on the target source. For the 345 GHz observations in 2007 and 
2009, after each 15 minutes we switched to 1733-130 (NRAO 



^ The Submillimeter Array is a joint project between the 
Smithsonian Astrophysical Observatory and the Academia Sinica 
Institute of Astronomy and Astrophysics and is funded b y the 
Smithsonian Institution and the Academia Sinica ( |Ho et al.|2004| ). 
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Table 1. Observational summary 



Date 


frequency 


array configuration 


baseline^ 


no. antennas 


integration time 


gain calibrator 




(GHz) 




(m) 




(min) 




12 Aug 2007 


345/355 


extended 


205 


7 


281 


1733-130 


10 May 2008 


681/691 +201/211 


compact north 


125 


5 


270 


Callisto 


6 Sep 2008 


681/691 + 345/355 


subcompact 


70 


6 


82 


Callisto 


7 Sep 2008 


681/691 + 345/355 


subcompact 


70 


7 


93 


Callisto, Ganymede 


4 July 2009 


345/355 


very extended 


516 


8 


330 


1733-130 



Notes. Maximum projected baseline 




19°5 



19°5 



19°5 



38^0 



345/355 GHz 



681/691 GHz 



-19°51'50" 



-19°51'55" 



38^0 



20 
15 
10 



Table 2. Beam sizes and noise levels 



o 
o 




150 



100 



o 
o 




o 
o 



Fig. 1. Continuum maps of GlO.47+0.03 observed with the 
SMA at 201/211 GHz (top panel), 345/355 GHz (central), and 
681/691 GHz (bottom). The map size is either 5 or W. Contours 
mark 20 and 50% of the peak flux, which is given in the lower 
right (in K). The beam is depicted in the lower left. The white 
crosses denote the Hii regions Bl, B2, and A (Cesaroni et al. 
|2010t from left to right). 



530), which lies 10.9° away from the target source. In addition, 
strong sources (Uranus, Ceres, 3c454.3, 3c273) were observed 
before and after these loops to allow bandpass and flux calibra- 
tion. 



Frequency 


beam"^ 


rms^ 




(GHz) 


r;°) 


(mJy/Beam) 


(K) 


199.9-201.8 


4.27 X 1.84; 57 


38 


0.15 


209.9-211.8 


4.14 X 1.95; 55 


37 


0.13 


201/211 cont. 


4.14 X 1.48; 60 


2.9 


0.014 


344.1-345.9 


0.45x0.29; 18 


55 


4.3 


354.1-355.9 


0.42 X 0.29; 22 


55 


4.4 


345/355 cont. 


0.42 X 0.28; 20 


12 


1.0 


680.3-682.2 


2.12x0.65; 37 


840 


1.6 


690.3-692.2 


2.08 X 0.60; 38 


820 


1.7 


681/691 cont. 


1.77x0.50; 42 


227 


0.67 



Notes. Different weighting schemes were used for continuum and 
line images. Major and minor axis and position angle (from north to 
east) are given. rms noise is given for a frequency width of 1 km s~^ 
in the line maps. It is evaluated from -10 to -5'' in both RA and Dec 
relative to the phase center. 



Calibration and editing of the data were done in the IDL 
MIR packag^ involving correction and application of the sys- 
tem temperature, phase-only bandpass calibration, continuum 
regeneration (spectral averaging), phase-and-amplitude band- 
pass calibration, flagging, gain calibration (phase and ampli- 
tude), and flux cal ibration. The data were converted to MIRIACj^ 
( Sault et al.|1995| ), where the edge channels of each chunk were 
flagged and the channels were rebinned to a width of 1 km s"^ 
at 201 GHz, 0.75 km s'^ at 345 GHz, and 1.42 km s'^ at 



681/691 GHz. Data from difl'erent dates were merged for each 
frequency setup. Since the 2007 observations had the phase cen- 
ter at R.A. 18:08:38.28, Dec. -19:51:50.0, while the later obser- 
vations were centered on R.A. 18:08:38.232, Dec. -19:51:50.4, 
the 345/355 GHz data set had to be merged in AIPS. The chan- 
nels less afl'ected by spectral lines were identified and used for 
separating lines and continuum in MIRIAD. The continuum data 
of upper and lower sideband were merged and imaged using 
almost uniform weighting at 201/211 GHz and 345/355 GHz 
(robust^ -2) and more natural weighting at 681/691 GHz 
(robust= 0.5). A cutoff' for cleaning of about 3 times the rms 
noise in the image was used. Based on these clean components, 
the continuum visibility data were self-calibrated, which sub- 
stantially reduced the ffuctuations outside the source, and the so- 
lutions were applied to the line data as well. Imaging of the spec- 
trally resolved data was done with robust= 0.5 at 201/211 GHz 
and 345/355 GHz and robust= 2 at 681/691 GHz. Table[2]gives 
the resulting beam sizes and noise levels in the maps. 



^ http ://w w w. cf a.harvard. edu/cqi/mircook.html 
^ http ://bima. astro . umd.edu/miriad 
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Line identification was made with the myXCLASS pro- 
gran]^ which accesses the and JPLpl molecular data 
bases. All figures in this paper were made with the GILDAS 
softwar^ 



3. Observational results 

3.1. Continuum 

Figure [T] shows the obtained continuum maps. The total flux 
is 6 Jy at 201/211 GHz, 27 Jy at 345/355 GHz, and 95 Jy at 
681/691 GHz, corresponding to a spectral index of 2.8 between 
the lower two and 1.8 between the upper two frequencies. While 
the beam sizes are not sufficient to resolve the continuum emis- 
sion at 201/211 and 681/691 GHz, the extension can be clearly 
seen at 345/355 GHz. At this frequency, the peak intensity is 158 
K (1.84 Jy/Beam), and the 3cr contour extends over y\ 

3.2. Line identification 

The data cubes were convolved to a common resolution of 5'' 
and central spectra were extracted. This resolution was chosen to 
decrease the noise and the eff'ects of a strong continuum, as well 
as to facilitate a comparison between the diff'erent frequencies; 
it gives the total (source-integrated) flux for most lines. 

To identify the spectral features, a simple homogeneous 
model was computed in Local Thermodynamic Equilibrium 
(LTE), using the myXCLASS program, which takes the difl'erent 
optical depths of the transitions into account. Synthetic spectra 
are compared to the data. A source size of 1.5'' (half-maximum 
diameter of a Gaussian source), a temperature of 200 K and a 
line width of 5 km s"^ were fixed, while the column density of 
each molecule was varied to obtain a good fit to the data (see 
Table [3|. The continuum was neglected in the radiative trans- 
fer, only a foreground absorbing column density of 5 x 10^"^ H2 
cm"^ was taken into account to approximate the higher absorp- 
tion at the higher frequencies (corresponding to r = 0.3 at 345 
GHz and r = 1.3 at 690 GHz). Although this is only foreground 
absorption, it serves to account for both absorption and contin- 
uum emission, as the latter weakens the lines as well. Continuum 
levels of 6, 8, and 7 K were added to the spectra at 201/211, 
345/355, and 681/691 GHz, respectively. Figure s [Al] [A^ and 
A.3 show the spectra and the model, and Table |3J gives the de- 
rived column densities. We note that the model is not supposed 
to give an optimum fit to the data, but is just for identification 
purposes. 



Table 3. List of identified molecules 



3.3. Line maps 

From many maps of interesting lines (see appendix. Figs. |B. If 
|B.5| ), we show here only a selection, which are relevant for the 
expansion motion (Fig.|2]) and the high excitation (Fig. [3]). They 
are ordered as in Table]3| The velocity was integrated over the 
ranges with detectable signal and no significant deviations in the 
channel maps - i.e. the whole flux for simple line shapes and sev- 
eral velocity ranges for more complex line shapes. For Fig.|2] we 
chose common velocity ranges for the four lines, which are 30- 
50 km s"\ representing the high- velocity part of the (front-side) 
outflow, 51-64 km s"^ the low- velocity part, 65-70 km s"\ the 





IIU. lillCIS 


Loiuiiiii ueii&iiy 






(cm 


SO 


5 


4(17) 


S02 


5 


3(17) 


ocs 


2 


3(18) 


CN 


1 


C 


HCN 


8 


1(18); 3(18)'' 


HNC 


1 


3(17) 


CH3NH2 


2 


4(17) 


CH3NC 


4 


7(15) 


NH2CN 


7 


2(16) 


NH2CHO 


15 


2(17) 


HC3N 


76 


5(16); 1(18)'' 


C2H3CN 


33 


7(17) 


C2H5CN 


55 


9(17) 


H2O 


1 


1(20) 


CO 


2 


C 


H2CO 


6 


3(18) 


CH3OH 


23 


9(18) 


H2C2O 


4 


4(17) 


C2H5OH 


23 


6(17) 


CH3OCH3 


14 


1.5(18) 


CH3OCHO 


38 


7(17) 


CH3CH3CO 


26 


5(17) 



Notes. The molecules are ordered by S-, N-, and 0-beari ng and b^ 



A.l 



^ https ://w w w. astro .uni-koeln. de/proj ects/schilke/XCLAS S 

^ http://www.cdms.de 

^ http://spec.jpl.nasa.gov 

^ http://www.iram.fr/IRAMFR/GILDAS 



complexity. Number of identified lines (labeled in Figs. 

Column density used in the myXCLASS model, with a source size 
of 1.5'' and a temperature of 200 K. The parentheses are powers of 10. 

The lines are in absorption, which is not modeled here. For HCN 
and HC3N, an additional component of 0.5'' and 500 K was used to 
more closely match the vibrational lines. 



systemic velocity, and 71-84 km s"\ the low- velocity part of the 
back- side outflow. Higher velocities are not detected, probably 
due to dust absorption. In case of contamination by a neighbor- 
ing line, the velocity range was chosen to avoid this line. Still, 
blending is likely for SO at 344.3 GHz by methanol (frequency 
corresponds to a 1.5 km s"^ lower velocity), H^^CN at 345.3 
GHz by SO2 (1 km s"^ higher velocity), H^^NC at 355.4 GHz by 
CH3CH3CO (3.3 km s-i lower velocity), and CH3NH2 at 354.8 
GHz by CH3OCHO (3.8 km s"^ higher velocity). 



3.4. Spectra 

Figure |4] shows central spectra of the 12 transitions with absorp- 
tion features. Most of the absorption is blue-shifted relative to 
the systemic velocity of 68 km s"^ 



3.5. Comparison to APEX data 

To estimate calibration and filtering of extended emission, we 
compare the SMA data to data from t he APEX (Atacama 
Pathfinder Experiment) 12-m telescope ( [Giisten et al. 20Q6[ 
ISchullere t al. 2009 ; Rolff's et al.|2011b|). The 345 GHz flux from 
the LABOCA bolometer array is 34 Jy/Beam, the 850 GHz flux 
from SAB OCA is 450 Jy/Beam. At both frequencies, the source 
is unresolved with the beams of 18.2 and 7.4'', respectively. The 
spectral index is 2.86, so the 690 GHz flux would be 250 Jy. The 
SMA fluxes at 210 and 345 GHz fit very well to the APEX data, 
but at 690 GHz the continuum is much lower. This might be 
partly explained by decorrelation due to fast phase fluctuations. 
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Fig. 2. Line maps (selected to trace expansion motions) integrated over the velocity range indicated in each panel (which can be 
considered, from left to right, as high- velocity part of the outflow, low-velocity part, systemic velocity, and low- velocity part of the 
outflow at the far side). The first velocity range of HC3N is blended by other lines, in the second velocity range of CO the emission 
is so extended to cause imaging artifacts. The map size is 5" (tick spaces are 1", centered on R.A. 18:08:38.236, Dec. -19:51:50.25). 
Beams are shown in the lower left, and the number in the lower right of each panel is the maximum flux in K km s~' (contours are 
±20 and 50% of that value). The color scale is from -1000 to 1000 K km s"' . The energy of the lower level is given in the upper left. 



■1,05^ 























-1 1 1 .l|207^ 


■ 



and partly by missing short baselines filtering out extended emis- 
sion. 

Figure [5] shows APEX spectra overlaid with SMA spectra 
which were convolved to the APEX beam (18.2'' at 345 GHz, 
17.7'' at 355 GHz, and 9.1" at 690 GHz). The lines match very 
well. The most notable exception is CO 6-5, which is purely in 
absorption with the SMA, but has strong emission as seen with 
APEX (>1000 Jy). Extended emission of CO 6-5 is filtered out 
due to missing short spacings, leaving an apparant absorption 
towards the continuum, which is more compact and hence not 
aff'ected as much as CO by the filtering. 



4. Modeling 

In this section, we compare a few models to continuum and vi- 
brationally excited HCN in order to constrain the structure of this 
complex region, mainly density distribution and velocity field. 
Testing diff'erent models of the source structure allows to exclude 
certain structures and to approach the real structure by adapting 
the models. The models are inevitably simplified, and we used a 
trial-and-error technique with parameter variations and compar- 
ison of model and data by eye to approach a good fit if possible 
for a tested structure, but no global optimization of fit parameters 
was performed. 
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H^^CN,V2=1 (345.2) 59-78 HCN,V2= 1 (354.5) 60-76 H^^CN,V2= 1 (690.3) 64-72 HC^^N,V2= 1 (69 1 .6) 63-73 




H^^CCCN,V7=2(354.8) 64-7 'HC3N,V4=V7= 1 (355. 1 ) 62-71 HC3N,V6= 1 (355.6) 64-80 H2^^0(692.1) 63-71 




Fig. 3. Selected high-excitation line maps. Above each panel, the molecule, the frequency of the transition (in GHz), and the velocity 
range (in km s"^) are given. The color scale ranges from to 1000 K km s"^ at 345/355 GHz, where the map size is 5'', and from 
to 100 K km s"^ at 690 GHz, where the map size is 10''. Beams are shown in the lower left, and the number in the lower right of 
each panel is the maximum flux in K km s"^ (contours are ±20 and 50% of that value). The energy of the lower level is given in the 
upper left. 



The three-dimensional radiative-transfer code RADMC- 
3E[^ developed by C. Dullemond, was employed to compute the 
dust temperature from stellar heating and the continuum and line 
emission of a model . The basic setup is the same as described in 



RolfFs et al. 



but here we also test dumpiness and ve- 
locity structure in the models. The line transfer assumes LTE, 
which is a good approximation for vibrationally excited HCN, 
but prevents the modeling of the ground- state lines. 



density fluctuations on small scales, which are expected from 
self-gravity and which improve the line fitting. Model E is an 
attempt to include the outflow (which is seen by the blue- shifted 
absorption features) and is presented in more detail. 

The models are described explicitly in the following. Table |4] 
summarizes the main properties of the models. 



4.1. The Models 

While no perfect fit was found, we selected five of the tested 
models for presentation and comparison to the data. Models A 
and B have been previou sly used to fit APEX a nd VLA data of 
this source, respectively ( [Rolffs et"ari|2011b|a| ). Model C has a 
radial density profile that best matches the observed continuum 
(if heated by the stars in the Hii regions). Models D and E have 



Model A This is the model that was presented in RolfFs et al. 
( 2011b| ) to fit the APEX data (continuum and many lines from 
HCN, HCO+, and CO, including vibrationally excited HCN). 



The density follows a radial power law, ^ = 7xl0^x( ) 



-1.75 



H2 cm"^ for radii larger than 485 AU. Inside there is an Hii re- 
gion with an electron density of 1.5 x 10^ cm"^ to reproduce 
the free-free radiati on from Bl, and a star of 5.6 x 10^ L©. The 
dust opacity is from'Qssenkopf & Henning| ( 1994| ) without grain 
mantles, but with coagulation at a density of 10^ cm"^. The HCN 

abundance is 3 x 10"^ at temperatures above 100 K and 5 x 10"^ 

http://www.ita.uni-heidelberg.de/~dullemond/software/radmc-3d below 100 K. The line width is 5 km s"^ FWHM, and the gas 
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Fig. 4. Lines with absorption features toward the continuum peak. The systemic velocity of 68 km s"^ and a blue-shifted velocity 
of 50 km s"^ are marked. There are absorption components centered at both velocities, indicating expansion. The dashed horizontal 
line denotes the continuum level. 



Table 4. Summary of the models described in Sect. 4.1 



Model 



density 



heating 



velocity field 



A power-law 1 star 

B Gaussian 3 stars 

C Plummer 3 stars 

D clumpy 3 stars 

E clumpy, outflow cone 3 stars 



infall 
none 
none 
dispersion 



displayed in Figs. 



_7 9|(red) 

^lm^j^(red) 
l^cyan). 



(green ) 

dispersion, outflow [7| (magenta), K ' 'll|14 (red) 



is radially infalling with 1 km s"^ to the center at Bl. The total 
mass (in a cube of 3 pc diameter) is 2.4 x 10^ M©, of which 70 
M© are at temperatures above 300 K. 



Model B This is the model that was presented in Rolflfs et al. 
P011a| ) to fit the VLA data of vibrationally excited HCN. The 
density follows a Gaussian, centered at B 1, with 7 x 10^ H2 cm"^ 
at the half-maximum radius of 7000 AU. As in all following 
models, heating sources are the stars in the Hii regions Bl with 
10^ L©, B2 with 8.3 x 10^ L©, and A with 4.6 x 10^ L©, which 
are placed in the plane of the sky. . The dust opacity is from 
[Ossenkopf & Henning| ( [1994| ) without grain mantles or coagu- 
lation, as in all following models. The HCN abundance is 10"^. 
The intrinsic line width (FWHM) is 8.3 km s"\ and there is no 
macroscopic velocity field. The total mass is 3.5 x 10^ M©, of 
which 350 M© are at temperatures above 300 K. 



Model C The density in this model follows a Plummer profile, 

n = 1.4x10^x^1-^ ( n5QQAu )^) ^2 cm"^ (half-maximum 
radius 6500 AU). It is centered 2000 AU south and 1000 AU 



west of B 1 . The Plummer model is very similar to a Gaussian 
inside the half-maximum radius, but has higher density outside 
(still falling oflf steeply as r"^). The HCN abundance is 10"^ at 
temperatures above 300 K and 10"^ below, as in the following 
models. The intrinsic line width (FWHM) is 8.3 km s"^ and 
there is no macroscopic velocity field. The total mass is 6.8 x 10^ 
M©, of which 350 M© are at temperatures above 300 K. 



Model D This model consists of 100 clumps, which have a 
Plummer half-maximum radius of 1000 AU and central densi- 
ties ranging from 4 x 10^ to 2 x 10^ H2 cm"^. The distribution of 
clump masses follows the same slope as the stellar Initial Mass 
Function. They are randomly placed in the model according to a 
Plummer distribution with half-maximum radius 4000 AU. The 
density in each model cell is the maximum of all contributions 
from the clumps. So with a lot of overlap, this structure consists 
rather of density fluctuations than of separate cores. The intrin- 
sic line width (FWHM) is 5 km s"^ Each clump has a random 
line-of-sight velocity with Gaussian half-maximum value of ±5 
km s"^ The line-of-sight velocity of each cell is an average of all 
clumps, weighted with their density contribution. The total mass 
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Fig. 5. Comparison of APEX data (green) to the SMA data convolved to the APEX beam (black). The dashed lines mark the 
continuum level as expected from LAB OCA and SAB OCA measurements. The HCN and CO lines and the 690 GHz continuum are 
affected by filtering of extended emission due to missing short spacings; the rest agrees very well. 



is 9.4 X 10^ Mq, of which 470 M© are at temperatures above 300 component, which cannot be exactly reproduced in diff'erent runs 



K. 



with the same parameters. 



Model E This model has the same structure as model D, but a 
half-maximum radius of the distribution of the clumps of only 
3000 AU. In addition, a bipolar outflow is present, with an open- 
ing angle of 100° and a length of 2 x 10^ AU. Its center Hes 1000 
AU south and 500 AU west of B 1 . The axis is tilted from the 
line of sight by 20°, and the foreground part directed towards 
B2. Inside the outflow cone, the density is reduced to 20% of 
the original value (thus preserving the clumpy structure) and the 
velocity is 10 km s"^ outwards. This is only a toy model of the 
outflow. The total mass is 4.6 x 10^ M©, of which 400 M© are at 
temperatures above 300 K. The clumpy structure is visualized in 
Fig. (6] 



4.2. Continuum 

The continuum radiation that the models emit is compared to 
the data. Models D and E have randomly placed core centers; 
therefore their radiation (continuum and lines) have a random 



The radial profile of the 345/355 GHz continuum map 
is extracted with central position R.A. 18:08:38.237, Dec. 
-19:51:50.421, which is the center of a two-dimensional 
Gaussian fit (1.4'' x 1.17'', elongated along the B1-B2 axis). 
The model maps were Fourier-transformed, folded with the uv- 
coverage of the observations, and imaged in the same way as the 
data. The radial profile of the model was extracted from the same 
central coordinates. Figure [7] shows a comparison of the radial 
profiles. Fig. [8] a comparison of the continuum map of model E 
to the data. Model A (power-law) fits worse to the radial profile, 
model C (Plummer) fits best. 

The flux from APEX/LABOCA (345 GHz) is 34 Jy/Beam 
(18.2" beam size) and from APEX/SABOCA (850 GHz) 450 
Jy/Beam (7.4" beam size). The models emit 31 (A), 20 (B), 26 
(C), 20 (D), and 18 (E) Jy/Beam for LABOCA and 266 (A), 231 
(B), 311 (C), 146 (D), and 162 (E) Jy/Beam for SABOCA. In 
models B-E an extended component similar to A could be added 
without afl'ecting the results of the interferometer modeling. 
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Fig. 11. Model E (right) compared to the /=13 direct ^-type line of vibrationally excited HCN at 40.8 GHz (VLA), shown as 
integrated line map and spectra at two locations (left, model spectra are overlaid in red). The white contours denote the 7 mm 
continuum in steps of 1000 K. 



4.3. Vibrationally excited HCN 

The models are compared to lines from vibrationally excited 
HCN, including the /=13 d irect ^-ty pe transition at 40.7669 
GHz observed with the VLA ( [Rolffs et a l. 201 la). Figures[9]and 
10 show this line and three rotational transitions of vibrationally 
excited HCN observed with the SMA, overlaid with models A- 
D. To reproduce the direct ^-type line, the rotational lines must 
be very optically thick, leading to self-absorption in the models. 
This is inevitable for optically thick lines from centrally heated 
spheres; the data however do not show self-absorption. The self- 



absorption decreases from model A to E. Model E is shown in 
more detail in comparison to the /=13 direct ^-type line of vi- 
brationally excited HCN (Fig. 11 ), the 4-3 line of vibrationally 
excited HCN (Fig.[l2]) and H^^N (Fig. [l3])^nd the 8-7 tran- 
sition of vibrationally excited H^^CN (Fig7[l4|). The outflow is 
imprinted in the line shapes. The levels of all these transitions lie 
between 1050 and 1400 K above ground. 
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SMA: HCN.V2=1 4-3 model for Gl 0.47+0.03 




Fig. 12. Model E (right) compared to the J=4-3 line of vibrationally excited HCN at 354.5 GHz (SMA), shown as integrated line 
map and spectra at two locations (left, model spectra are overlaid in red). The white contours denote the 355 GHz continuum in 
steps of 25 K. 



SMA: h"cn,V2=1 4-3 model for Gl 0.47+0,03 




Fig. 13. Model E (right) compared to the /=4-3 line of vibrationally excited H^^CN at 345.2 GHz (SMA), shown as integrated line 
map and spectra at two locations (left, model spectra are overlaid in red). The white contours denote the 345 GHz continuum in 
steps of 25 K. 
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Fig. 14. Model E (right) compared to the /=8-7 line of vibrationally excited H^^CN at 690.4 GHz (SMA), shown as integrated line 
map and spectra at two locations (left, model spectra are overlaid in red). The white contours denote the 690 GHz continuum in 
steps of 25 K. 




Fig. 6. Clumpy structure of model E displayed as isocontours of 
10^ H2 cm"^ (green, transparent). The Hii regions, whose stars 
heat up the gas, are shown in blue. The size of the contoured 
region is around 0.2 pc. 

5. Discussion 

5.1. Chemistry 

In addition to about 350 identified lines, there are around 90 
unidentified (U-) lines. 75 of these 90 U-lines are in the 201/211 
GHz range, which covers a similar velocity range as the higher- 



T — I — I — r-i—| 1 1 1 1 — I — I — r-i—| 1 r 




Radius [arcsec] 

Fig. 7. Models A (red), B (green), C (blue), D (cyan), and E (ma- 
genta) compared to the SMA continuum profile, which is dis- 
played as data points with a O.T' radial interval. The errorbars 
denote the rms deviation from a circular shape, and the dashed 
curve represents the beam. 



frequency bands combined (almost 6000 km s"^). This is proba- 
bly due to the lower noise at these frequencies. The line density 
is mainly determined by the width of single lines (5-10 km s"^), 
i.e. a forest of overlapping lines covers most of the bands. The 
identified molecules (3 S -bearing, 10 N-bearing, and 9 other. 
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Fig. 9. Models A (red) and B (green) compared to lines from Fig. 10. Models C (red) and D (green) compared to lines from 
vibrationally excited HCN at the position of B 1 and 0.3'' (VLA) vibrationally excited HCN at the position of B 1 and 0.3'' (VLA) 
or 0.8" (SMA) south of B L or 0.8" (SMA) south of B L 



0-bearing, molecules) represent only the strongest lines, while 
more complex molecules have many weak lines which probably 
add up to the continuum and are undetectable due to line confu- 
sion. No spatial separation of different molecules could be seen. 
The differences in the maps can probably be explained by exci- 
tation and optical depth. 

In the detailed models (Sect. |4]), the abundance of HCN in 
the hot gas is very high, on the order of 10"^ relative to H2. This 
is needed to reproduce the direct ^-type line and the vibrationally 
excited H^^CN. Around 5% of the nitrogen and 1.4% of the car- 
bon is in HCN (for solar N/H of 1.1 x lO'^ and C/H of 3.6 x lO'^, 
[Anders & G revesse 1989 ), and HCN/CO is on the order of 0.1 
in the hot central region. We note that such a high HCN abun- 
dance needs an explanation by future chemical models, incor- 
porating high-temperature reaction networks. The abundance of 
HC3N seems to follow a similar increase at high temperatur es, a s 
the s trong emission from vibrational states suggests (Figs.|A.l| 



which at high temperatures re acts to form water instead of de- 
stroying N-bearing molecules ( [Rodgers & Charnley|2001 1. 



A3]). Notable line detections are HC^^N,V2=1 at 691.6 GHz 
(lower level is 1139 K above ground), HC3N,V4=V7=1 at 355.1 
GHz (1908 K), and the ^^C substituted HC3N,V7=2 (970 K). The 
H^^O line at 692.1 GHz (694 K) implies a high abundance of 
water in the dense, warm gas. These abundance enhancements 
could be connected, e.g. through a reduced abundance of OH, 



5.2. Density Distribution 

As the high-resolution continuum map at 345/355 GHz shows, 
a power-law radial density is not consistent with the data. A 
Gaussian fits well in the inner part, but falls off too steeply in 
the outer parts. What fits best to the data is a Plummer model 
(Fig. [7] blue). This profile is very similar to a Gaussian inside 
the half-maximum radius, but is denser outside, where it falls off 
as r"^ . It is also used to describe the stellar density of star clus- 
ters. Although that might reveal a connection of this forming star 
cluster to later stages, it could also have completely different 
physical reasons. The stellar density distribution is determined 
by stellar dynamics, while the gas density distribution is deter- 
mined by the interplay of different sources of pressure (gravita- 
tional versus turbulent, rotational, magnetic, radiative, and ther- 
mal pressure). The observed central flattening of the density can 
be explained by centrally increased pressure, which is expected 
from feedback by the newly formed massive stars. This pressure 
stalls the infall and piles up the mass where the infall stops, re- 
sulting in a less centrally peaked mass distribution. 
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The SMA continuum at 690 GHz is rather weak compared 
to the 345 GHz flux and the high-frequency single-dish flux. As 
interferometer observations at these frequencies a re challenging 
and only very few have been conducted yet (e.g. [Beuther et al.| 
|20Q6| ), this could be due to decorrelation caused by rapid phase 
noise. On the other hand, the line strengths are similar to the 
APEX data (Fig. [5]), supporting the data quality and the cali- 
bration, so if the low continuum is real, it can be due to opti- 
cal depth efl'ects in very dense, small condensations, while ex- 
tended emission is filtered out. The continuum at 690 GHz must 
be more extended than at 345 GHz because the larger dust opac- 
ity makes it stronger only in the optically thin regions, while it 
can be weaker in the optically thick case - instead of the whole 
hot interior, one sees only up to the colder foreground. In addi- 
tion, the image fidelity at 690 GHz is not as good as at 345 GHz, 
in particular the shortest baselines are twice as large. As the flux 
rises with shorter baselines, this means much more filtering out 
of extended emission at 690 GHz. 

There is no fragmentation observed in neither the continuum 
nor the lines. The reason may be that we do not resolve struc- 
tures smaller than the beam size of about 3000 AU. We note that 
the continuum source extends over about 10 beams (3''), so no 
large-scale fragmentation is present. Models of the line emission 
reveal that a spherical density distribution leads to strong self- 
absorption features, which are not observed. Clumpiness better 
reproduces the line shapes, but without sufficient spatial resolu- 
tion we cannot derive the properties of fragmentation. 

This missing fragmentation is sim ilar to the contiuum data 
recently obtained by Qin et al. ( 2011| ) for SgrB2-N, which fol- 
lows a very similar radial profile and also displays a similar spec- 
tral line content. Maybe the two sources are at the same evo- 
lutionary stage, where the first hypercompact Hn regions have 
developed and heated the gas, but large-scale fragmentation has 
not set in. 



The rotational transition of vibrationally excited HCN (4-3 
at 354.46 GHz) has a very high optical depth. The Einstein A 
coefficient is 50,000 times larger than for the direct ^-type line 
(/=13 at 40.7669 GHz, observed with the VLA), which has an 
optical depth of < 1. Also the comparison to the H^^CN line 
(Fig. [3]) shows its high optical depth. In models which are ho- 
mogeneous on small scales (smooth density gradient, such as 
models A-C), this optically thick line is self-absorbed due to 
the temperature gradients caused by internal heating. Such large 
temperature gradients are inevitable as the stars must be deeply 
embedded, shown by the high efficiency of heating and the com- 
pactness of the Hii regions. The observed single-peaked profile 
can only be explained by an inhomogeneous density and velocity 
field (such as models D-E), where the self- absorption is smeared 
out. Such small-scale clumpiness might persist in later stages of 
evolution, when the gas is ionized ( Ignace & Churchwell|2004 l. 



The observed continuum depends on density, temperature, 
and dust opacity. The density can vary by many orders of mag- 
nitude, and therefore dominates the resulting continuum. It is 
clear that the temperature is high in the inner region and falls off' 
outwards; this is included in the models as the stars in the hyper- 
compact Hn regions heat the dust. It is entirely possible, how- 
ever, that more heating sources are present, although they are 
not needed for heating the central region. The extension of high- 
excitation lines is generally a bit smaller in the models than in 
the data, especially highly excited HC3N is somewhat extended 
to the east (Fig.[3|. 



5.3. Velocity Field 

While most lines do not show any velocity structure in their 
maps, there seems to be a component of highly excited CH3OH 
at lower velocity close to the UCHn region A (Fig. [3]). 
Alternatively, this could arise from blending by a diff'erent line. 

On large scales, the asymmetries of self-absorbed lines 
clearly indicate infall motions (Rolff's et al. 2011b, and Fig. [5]), 
and also the absorption of the H2CO line is a bit red-shifted 
(Fig. |4]). On smaller scales, we see expansion motions, as very 
nicely traced by the absorption features (Fig. [4]). The blue-shifted 
absorption is stronger in the western part of the continuum. In 
addition, the SO line at 344.3 GHz, the HCN fine at 354.5 GHz, 
and the CO line at 345.8 GHz have red-shifted emission in the 
north-eastern part. A red-shifted clump 5'' north of the hot core, 
which is only seen in HCN and CO, seems to be too far away to 
be related to the outflow (Fig. [2]). 

l Olmi et an([T9"96) detected a north-so uth velocity gradient in 



Hofner & Churchwell (1996]) see blue- shifted water 
north and 0.4'' south of Bl, and a red- shifted water 



i^CO, and 
masers 0.2 

maser (at 70 km s"^) 1.2'' north of Bl. Absorption against the 
hypercompact Hii regions Bl and B2 is found to be blue-shifted 
in NH3(4,4) (Cesaroni et al.,2010, 53 km s"^ toward Bl and 61 
km s"^ toward B2), and also partly in vibrationally excited HCN 
( [Rolff'setal.|2011a l. and Fig.[TT]). 

One can conclude that there is likely a bipolar outflow, whose 
foreground part is tilted to the south-west, and whose back- 
ground part is tilted to the north-east. A sketch of the outflow 
scenario is presented in Fig. 15 This scenario explains qualita- 
tively all observed velocity features. The difl'erent NH3(4,4) ve- 
locities, for instance, can be understood as an efl'ect of the pro- 
jection of the outflow to the line-of-sight. An embedded source, 
which has not yet developed a detectable Hii region, is driving 
the outflow in this scenario. A nearly spherical infall with lower 
velocity, but higher rate than the expansion surrounds the out- 
flow and produces the blue asymmetric line profiles seen on 
larger scales. This accretion probably continues well into the 
central region along certain paths, driven by gravitational attrac- 
tion and the momentum of the large-scale accretion ffow. 

An alternative scenario might arise if the hypercompact Hn 
regions Bl and B2 are not spherical, but in the form of bipolar 
expanding bubbles from a photoevaporating disk ( Hollenbach] 
etar||T9^ . In this case, the stars in the HCHn regions would 



drive two outflows, which are both roughly aligned with the line- 
of-sight. 

In addition, it is possibile that the source is a large pseudo- 
disk accreting from a massive envelope, rotating roughly in the 
plane of the sky, driving an outflow perpendicular to it, and hav- 
ing formed Bl and B2 as binaries (and maybe more stars). This 
would correspond to the transition between phase II and III of 
Zapata et al.] ( |2010 ). 



A quantitative radiative transfer model, which reproduces all 
observed lines, is possible, but beyond the scope of this paper. 
That the outflow is important for reproducing the line shapes 
has been demonstrated in our modeling, as only model E could 
reproduce the shapes (Sect.|4]). 



6. Conclusions 

With the SMA, we have obtained high-resolution, spectrally re- 
solved maps of the massive hot molecular core G10.47-h0.03 at 
difl'erent frequencies, covering a bandwidth of 12 GHz in total 
and including observations at 690 GHz. Our main results are 
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observer 




Fig. 15. Qualitative scenario for the outflow in GlO.47+0.03. 
The observer is above the plot. Infalling molecular gas is shown 
in green, expanding molecular gas in blue, and ionized gas in 
white. The circles denote half-maximum and Plummer radius in 
Model C (the gas must be clumpy, though) and the boundaries 
of the Hii regions, which have density gradients ( [Cesaroni et al. 
2010j). If the Hii regions are not spherical, the outflow could also 
originate from the ionizing stars. 
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- Hundreds of molecular lines reveal a rich chemistry, with 
molecules such as HCN and HC3N especially abundant at 
high temperatures. This high abundance demands explana- 
tions from chemical models. Vibrationally excited HCN, 
whose levels lie more than 1000 K above ground, is very 
optically thick, and even vibrationally excited HC^^N shows 
up at 690 GHz. 

- Blue-shifted absorption in a dozen lines indicates an outflow 
oriented roughly along the line-of- sight. It is embedded as 
there is also absorption at the systemic velocity of cold fore- 
ground gas in CO and HCN. 

- The averaged radial profile of the submm continuum dis- 
plays a central flattening and rapid fallofl', and is best fitted 
by a Plummer profile of the density. The mass of the source 
is on the order of several thousand M©, of which a few hun- 
dred are at high (>300 K) temperatures. 

- No fragmentation is observed over the ~10 beam sizes 
(30,000 AU) of the continuum emission. High-excitation 
lines which are very optically thick do not show self- 
absorption. The line modeling suggests that this must be due 
to density fluctuations in combination with the velocity field. 
If the rather low continuum at 690 GHz is real, it points to 
small, very dense condensations. 

From these findings, a picture emerges of a very young form- 
ing star cluster, characterized by beginning of feedback from 
massive stars, while infall is ongoing. This feedback includes 
heating of the dust and gas, ionization that is still confined to 
small regions, and increased pressure in the inner part that leads 
to expansion motions and central flattening of the density. 
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Appendix A: Full Spectral Range 



In Figs. |Al] |A3| and [A3 

range, convolved to a reso! 



we show the full observed spectral 
ution of 5'', and label the identified 



lines with the molecule (see also Sect. |3.2| 



Appendix B: Line maps 

In Figs.|B]T]|B]2]|B31|B4| and|B3]we show 100 integrated line 
maps, ordered by molecule as in Table [3] and by frequency of the 
transition (see also Sect.|3.3|). 
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Fig. A.l. Full spectral range towards GlO.47+0.03 convolved to 5'' resolution. The green overlay is the myXLASS model for all 
molecules except HC3N, whose emission is overlaid in red. The displayed panels cover 199.85-201.85 GHz (upper three) and 



209.85-210.55 GHz (bottom); the rest is shown in Figs. A.2 and A.3 
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Fig.A.2. As Fig. 



A.l 



displaying 210.5-211.85 GHz (upper two) and 344.05-345.95 GHz (lower two). 
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Fig. A.3. As Fig. |AJj displaying 354.1-355.95 GHz (upper two), 680.2-682.3 GHz (third panel from top), and 690.2-692.3 GHz 
(bottom). 
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Fig.B.l. Integrated line maps. Above each panel, the molecule, the frequency of the transition (in GHz), and the velocity range (in 
km s"^) are given. The map size is either 5, 10, or 15'', as written in the map (tick spaces are V\ centered on R.A. 18:08:38.236, 
Dec. -19:51:50.25). Beams are shown in the lower left, and the number in the lower right of each panel is the maximum flux in 
K km s"^ (the contours are ±20 and 50% of that value). The color scale ranges from the minimum to the maximum value. The 
energy of the lower level is given in the upper left. 
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Fig.B.2. As Fig. B.l contd. 
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Fig.B.3. As Fig. B.l contd. 
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Fig.B.4. As Fig. B.l contd. 
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Fig.B.5. As Fig. B.l contd. 



